CMB constraints on Dark Matter models with large annihilation cross-section 
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The injection of secondary particles produced by Dark Matter (DM) annihilation around red- 
shift ~ 1000 would inevitably affect the process of recombination, leaving an imprint on Cosmic 
Microwave Background (CMB) anisotropies and polarization. We show that the most recent CMB 
measurements provided by the WMAP satellite mission place interesting constraints on DM self- 
annihilation rates, especially for models that exhibit a large Sommerfeld enhancement of the anni- 
hilation cross-section, as recently proposed to fit the PAMELA and ATIC results. Furthermore, we 
argue that upcoming CMB experiments such as Planck, will improve the constraints by at least one 
order of magnitude, thus providing a sensitive probe of the properties of DM particles. 
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I. INTRODUCTION 

The recent measurements of the Cosmic Mi- 
crowave Background (CMB) anisotropy and polariza- 
tion from experiments as WMAP ACBAR and 
BOOMERANG |[ have confirmed the theoretical pre- 
dictions of the standard cosmological model based on in- 
flation, dark matter and a cosmological constant. This 
not only permits to place stringent constraints on several 
parameters of the model but also to use it as new labo- 
ratory where to test physical processes in a environment 
not achievable otherwise. 

In particular, there is a remarkable agreement between 
the theoretical description of the recombination process, 
occurring at z r ~ 1000, and CMB data, a circumstance 
that severely constrains new sources of ionizing photons, 
and more in general any deviation from standard recom- 
bination [H, as recently shown by several groups of au- 
thors (see e.g. [||, @, Q). Most of the recent literature 
has analyzed the modified recombination by means of a 
phenomenological approach, parameterizing in a model 
independent way the modifications to the recombination 
process. Here, we focus instead on the constraints that 
can be placed on the properties of DM particles, under 
the assumption that recombination is modified only by 
dark matter annihilation. With respect previous studies 
@> our analysis includes more recent data (WMAP 
5-year data), and it concentrates on a new class of DM 
models that have been recently proposed to explain the 
observed anomalies in cosmic ray data. 

In fact, the attempt to explain the high energy positron 
and electron rise seen by PAMELA [l(| and ATIC [U in 
terms of Dark Matter (DM) annihilation has prompted 
the proliferation of new DM candidates with very large 
annihilation cross-section. In particular, in models with 
a "Sommerfeld" enhancement of the annihilation cross- 
section (av) , the efheent exchange of force carriers at low 



relative particle velocities leads to a velocity-dependent 
(av), which behaves roughly as cx 1/v for high v, and 
saturates below a critical v s (typically smaller than the 
local velocity dispersion, vq, see below), that depends on 
the ratio between the masses of the force carrier and the 
DM particle. A nice feature of these models is that they 
can be made naturally consistent with standard ther- 
mal freeze-out. In fact, DM freezes out typically with 
f3 = v/c = 0(1), and {crv) will grow from this minimum 
value as the universe cools and expand. Then, when the 
first gravitationally bound structures form, DM virializes 
within the gravitational potential of the host halo, thus 
leading for Milky Way (MW)-likc galaxies at z=0 to viri- 
alized velocities of order (3 ~ 10 . Smaller velocities, 
thus a larger (av) can be achieved in DM haloes with 
low velocities such as MW subhaloes[l2j or small haloes 
at high redshift 13]. although it is unclear whether the 
annihilation flux can be boosted enough to explain the 
PAMELA and ATIC data without being in conflict with 
other measurements, such as the anti-proton or gamma- 
ray fluxes towards the Galactic center (see e.g. [1J| and 
references therein). 

When recombination occurs, around z r ~1000, the 
relic WIMPs have not yet formed sizable gravitationally 
bound structures and are cold enough for the Sommer- 
feld mechanism to produce substantial enhancement of 
the annihilation cross-section with respect to the thermal 
value (after kinetic decoupling DM particle temperature 
evolves adiabatically as Toe z 2 , so 0(z r ) ~ 10~ 8 , for a 
O(100GeV/c 2 ) mass WIMP). As we will see below, the 
actual enhancement is model-dependent, because differ- 
ent DM models lead to a different behaviour of (av) z , but 
in general we expect that for large enough cross-sections, 
DM annihilation will significantly modify the recombina- 
tion history, thus leaving a clear imprint on the angular 
power spectra of CMB anisotropy and polarization. 

Our paper is organized as follows: in the next section 
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we describe the effects of annihilating dark matter on 
primordial recombination and the characteristic imprints 
on the CMB angular spectra. In section III we describe 
our analysis method. In section IV we present the results 
of our analysis. Finally, in Section V, we discuss our 
conclusions. 



II. ANNIHILATING DARK MATTER AND 
THERMAL HISTORY OF THE UNIVERSE 

Annihilating particles affect the ionization hystory of 
the Universe in three main different ways. The interac- 
tion of the shower produced by the annihilation with the 
thermal gas can i: ionize it, ii: induce Ly-a excitation 
of the hydrogen and Hi: heat the plasma; the first two 
modify the evolution of the free electron fraction x e , the 
third affects the temperature of baryons. In the presence 
of annihilating particles, the evolution of the ionization 
fraction x e satisfies: 



dx e 
dz 



1 



(l + z)H(z) 



[R s (z)-I s (z)-I x (z)], (1) 



where R s is the standard recombination rate, I s the ion- 
ization rate by standard sources, and Ix the ionization 
rate due to particle annihilation. 

Following the seminal papers [TEl ]. standard recombi- 
nation is described by: 

[R 8 (z) - I s (z)] = Cx [x 2 e n u a B - - x c )e- h ^' k * T *} 

where nu is the number density of hydrogen nuclei, 
qb and /3b are the effective recombination and photo- 
ionization rates for principle quantum numbers > 2 in 
Case B recombination, V2s is the frequency of the 2s level 
from the ground state and Tb is the temperature of the 
baryonic gas and the factor C is given by: 



C = 



[l + M 2s i s n H (l - x e )] 



[1 + KK 2sls n n {l - x e ) + Kf3 B n H (l - x c )] 



(3) 



where Ai S 2s is the decay rate of the metastable 2s level, 
nn(l— x e ) is the number of neutral ground state H atoms, 
and K = X^/(8ttH(z)) with H{z) the Hubble expansion 
factor at rcdshift z and A Q is the wavelength of the Ly-a 
transition from the 2p level to the Is level. 

The Ix term of equation [1] represents the contribute 
to the electron fraction evolution by a "non-standard" 
source; in our case it takes into account that during re- 
combination particle annihilation increases the ionization 
rate both by direct ionization from the ground state, and 
by contributing additional Lyman-a photons. The latter 
boosts the population at n = 2, increasing the rate of 
photoionization by the CMB from these excited states. 
Therefore, the ionization rate due to particle annihilation 



where Ixi is the ionization rate due to ionizing photons, 
and Ixa the ionization rate due to additional Lyman 
alpha photons. 

The rate of energy release per unit volume by a 
relic self-annihilating dark matter particle is given by 



(5) 



with tidm{z) being the relic DM abundance at a given 
redshift z, (o~v) is the effective self-annihilation rate and 
m x the mass of our dark matter particle, £Idm is the dark 
matter density parameter and p c the critical density of 
the Universe today; the parameter / indicates the frac- 
tion of energy which is absorbed overall by the gas, under 
the approximation the energy absorption takes place lo- 
cally. This on-the-spot approximation has been adopted 
by previous analysis (@) 

Each of the terms in Eq. Q] is related to the rate of 
energy release as: 



Ixi = C 

I Xa = (1-C)X, 



[dE/dt] 
n,H(z)Ei 
[dE/dt] 
' nH{z)E a 



(6) 

(7) 
(8) 

where Ei is the average ionization energy per baryon, E a 
is the difference in binding energy between the Is and 
2p energy levels of a hydrogen atom, n# is the number 
density of Hydrogen Nuclei and \i — Xa = {l — x e )/3 are 
the fractions of energy going to ionization and to Lyman- 
alpha photons respectively, given by I16H following the 
approach of Shull and Van Steenberg fl71 ]. 

Finally, a fraction of the energy released by annihilat- 
ing particles goes into heating of baryonic gas, adding an 
extra Kh term in the standard evolution equation for the 
matter temperature T&: 

m = 8a T a R T* x _ 

dz 6m e cH{z) 1 + /ho + x e 

+ 2T fc , (9) 



Ix(z) = Ixt(z) + Ixa(z), 



(4) 



ik B H{z) l + fn c + x 
where the non standard term is given by: 

(dE/dt) 

K h = Xh r-r- (10) 

n H (z) 

and Xh = (1 + 2x e )/3 is the fraction of energy going into 
heat given by [l6| . 



III. ANNIHILATING DARK MATTER AND 
THE CMB 

We can now compute the theoretical angular power in 
presence of DM annihilations, by modif ying the REC- 
FAST routine ([18]) in the CAMB code (|ljf), following 
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FIG. 1: Evolution of the free electron fraction as function of 
redshift for different values of p ann = [0, 1(T 6 , 5 x 10" 6 , 10~ 5 ] 
m 3 /s/Kg. 
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FIG. 2: TT, TE, EE angular power spectra (from Top to 
Bottom) for different values of p ann = [0, 1CT 6 , 5 x 1(T 6 , 1CT 5 ] 
m 3 /s/Kg . 



the prescription described in the previous section. The 
dependence on the properties of the DM particles is en- 
coded in the quantity 



< (TV > _ 
J — Pann 



(11) 



appearing in eq. [5l that we use as a parameter in the 
code. 

In Fig[T]we show the evolution of the free electron frac- 
tion for different values of p ann - As we can see, the DM 
annihilation model we consider can increase the free elec- 
tron fraction after z ~ 1000 by one order of magnitude, 
increasing the optical depth to last scattering surface and 
smearing the visibility function. The consequences of 
such annihilation can be seen in Fig[2]where we show the 
CMB anisotropy, cross-polarization and polarization an- 
gular power spectra for different values of p an n- DM an- 
nihilation damps the acoustic oscillations in the angular 
power spectra as in the case of an instantaneous reioniza- 
tion. However, large scale polarization is left unchanged 
by dark matter annihilation and a degeneracy between 
these two effects can indeed be broken. Although DM 
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FIG. 3: Constraints on the u>b, n s and u c parameters in the 
case of standard recombination (solid line), or including dark 
matter annihilation (dashed line). 



annihilation could play a role in the subsequent reion- 
ization of the Universe , the effect is likely to be small 
[2fj| |. unless one invokes very high anihilation cross sec- 
tions [2l|. Here, we don't consider a particular model for 
reionization, and simply adopt the parametrization of a 
full and instantaneous reionization at redshift z r < 30. 

We search for an imprint of self-annihilating dark mat- 
ter in current CMB angular spectra by making use of the 
publicly available Markov Chain Monte Carlo package 
cosmomc [22| . Other than p ann we sample the following 
six-dimensional set of cosmological parameters, adopting 
flat priors on them: the physical baryon and CDM densi- 
ties, uj\y = fl h 2 and u> c — fl c h 2 , the scalar spectral index, 
n s , the normalization, \nl0 10 A s (k — 0.05/Mpc), the op- 
tical depth to reionization, r, and the ratio of the sound 
horizon to the angular diameter distance, 9. 

We consider purely adiabatic initial conditions. The 
MCMC convergence diagnostic tests are performed on 4 
chains using the Gelman and Rubin "variance of chain 
mean" / "mean of chain variances" R— 1 statistic for each 
parameter. Our 1 — D and 2 — D constraints are ob- 
tained after marginalization over the remaining "nui- 
sance" parameters, again using the programs included 



Experiment p a nn 95% c.l. 



WMAP 


< 2.0 x 10 


- b m a /s/kg 


Planck 


< 1.5 x 10 


-'mVs/kg 


CV1 


< 5.0 x 10" 


s m^/s/kg 



TABLE I: Upper limit on p an n from current WMAP obser- 
vations and future upper limits achievable from the Planck 
satellite mission and from a cosmic variance limited experi- 
ment. 
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FIG. 4: Constraints on the self-annihilation cross-section at 
recombination (av) Zr times the gas-shower coupling param- 
eter /. The dark blue area is already excluded by WMAP5 
data, whereas the more stringent limit (dashed area) refers to 
the constraints which will be possible to apply with Planck. 
The light blue area is the zone ultimately allowed to probe by 
a cosmic variance limited experiment with angular resolution 
comparable to Planck. 



in the cosmomc package. We use a cosmic age top-hat 
prior as 10 Gyr < to < 20 Gyr. We include the five-year 
WMAP data t lJ (temperature and polarization) with 
the routine for computing the likelihood supplied by the 
WMAP team (we will refer to this analysis as WMAP5). 



IV. RESULTS 

Using the WMAP-5 dataset and applying the analysis 
method described in the previous section, we found an 
upper limit p ann < 2.0 x 10~ 6 m 3 /s/kg at 95% c.L, with 
no indications for modified recombination in agreement 
with previous and similar analyses. The implications of 
this limit are discussed in the next section. While we de- 
tect only an upper limit it is interesting, from a cosmo- 
logical point of view, to investigate the possible impact 
of this parameter on the estimation of other parameters 
as the baryon density Wb, the cold dark matter density 
uj c and the scalar spectral index ns- In Figure [3] we plot 
the 1-D likelihood distributions for these three parame- 
ters derived assuming the standard case (i.e. p ann = 0) 
and letting this parameter to vary freely. As we can see, 
including p ann into the analysis changes the constraints 
of uj b = 0.0228 ± 0.0006 and n s = 0.965 ± 0.014 (ob- 
tained in the standard case) to lo^ — 0.0230 ± 0.0006 and 
n s = 0.977 ± 0.018. The cosmological constraints on the 
cold dark matter density are on the contrary not affected 
by the inclusion of p a nn ■ 

With the advent of the Planck satellite mission, it is in- 
teresting to forecast to what extent the above limit will be 
improved by this mission. We have therefore forecasted 



future constraints on p ann assuming simulated Planck 
mock data with a fiducial model given by the best fit 
WMAP5 model (with standard recombination) and ex- 
perimental noise described by (see [HI): 



w 



-1/2 



/iK-rad 



exp 



FWHM 



/rad) 2 



Sin 2 



(12) 



with w 1 / 2 — 63/i-fT as the temperature noise level (we 
consider a factor \J~2 larger for polarization noise) and 
#fwhm = 7' for the beam size. We take f s k y = 0.65 
as sky coverage. We found that the Planck mission in 
the configuration described above will have the ability of 
placing a constraint of p an n < 1-5 x 10~ 7 m 3 /s/kg at 95% 
c.L 

It is also interesting to investigate the ultimate ability 
of cosmology to place constraints on p ann . We have there- 
fore repeated the analysis with an ideal Cosmic Variance 
Limited experiment with resolution up to £ ma x — 2500. 
In this case we found p a „„ < 5.0 x 10~ 8 m 3 /s/kg at 95% 
c.L 

These constraints are summarized in fig. 21 where we 
show the allowed values of f(crv) as a function of the 
WIMP mass m x , for the different experiments described 
above. These results place useful constraints on the DM 
annihilation cross-section at very small relative velocity. 
This is particularly important for models with a large 
"Sommerfeld enhancement" (SE), a non-perturbative ef- 
fect arising from the distortion of the wave functions 
of the two annihilating particles, due to the exchange 
of Coulomb-like forces mediated by (possibly new) force 
carriers [24| • The interest in these models arises from the 
fact that larger-than-thermal annihilation cross-section 
are required if one wants to explain the the rise in the 
electron and positron spectra observed by PAMELA and 
ATIC in terms of DM annihilation (see e.g. the discus- 
sion in Ref. [Hj]). We briefly recall here the basics of the 
SE. For two DM particles undergoing s-wave annihila- 
tion, the wave function in the non-relativistic limit obeys 
the Schrodingcr equation 

Tp"{r) - m x V{r)%P{r) + m^p 2 ip(r) = (13) 

In the limit where the mass of the carrier and the relative 
velocity of DM particles are small, it is easy to find an 
analytic approximation to the SE 



exp 



-air/{3^ 



(14) 



which exhibits the S ~ 1//3 behaviour that we mentioned 
in the introduction. Interestingly, a full calculation shows 
that the true solution saturates at (3 ~ m ( f > /m x , and it 
actually develops resonances, that lead to very large SE 
for specific combinations of masses and m x , and the 
coupling a. In order to compare the constraints on p ann 
obtained from the analysis of CMB data with theoreti- 
cal models, we have numerically integrated equation 1131 
assuming a Yukawa potential V(r) — — exp[— m,pr]a/r 
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FIG. 5: Constraints on the self-annihilation cross-section at recombination (av) Zr , assuming the gas-shower coupling parameter 
/=0.5, see text for details. Regions above the solid (/long dashed/short dashed) thick lines are ruled out by WMAP5 (/Planck 
forecast /Cosmic Variance limited); the thin dotted and dashed-dotted lines are the predictions of the "Sommerfeld" enhanced 
self-annihilation cross sections with force carrying bosons of ni0 = lGeV/c 2 and ni0=9OGeV/c 2 respectively, see text for details. 
Notice that these constraints apply to (av) at very low temperatures such that it is in saturated Sommerfeld regime, and 
therefore directly comparable with results from galatic substructures and dwarf galaxies constraints as from 



mediated by a boson of mass = lGeV/c 2 , and 
= 90GeV/c 2 , taking tv=l/47r (see e.g. [26| for 
details); /3=10 -8 as appropriate for the recombination 
epoch. 

The results are visualized in fig 03 were we show the 
region excluded by our analysis in the (av) vs. m x plane, 
corresponding to the 95 % c.l. upper limit on the cross 
section that can be derived by combining cq. [TT]with the 
constrains on p ann in table [TJ 



u z r .26 



= 71.2 



„max 
Pann 



(- 



WCL5 



^.O-lO-^s-Mcg-V VIOOGeV/ V / 

. (15) 

where avf 1 ^ denotes the upper limit of the annihilation 
cross section at recombination in units of 10 - cm s _1 . 

We have adopted in this formula, and in fig0 a fiducial 
value / = 0.5 for the coupling between the annihilation 
products and the gas, following the detailed calculation 
of DM-induced shower propagation and energy release 
performed by [27l |; this number is a good approximation 
averaged on all channels, although its actual value will 
ultimately depend on the composition of the shower and 
on its energy spectrum, namely on the nature of the anni- 
hilating DM particle itself. It is however straightforward 
to obtain the constraints for different values of /. 



We find that the most extreme enhancements are al- 
ready ruled out by existing CMB data, while enhance- 
ments of order 10 3 -10 4 with respect to thermal value 
(av)=3x 10~ 26 cm 3 /s, required to explain the PAMELA 
and ATIC data, will be probed over a larger WIMP mass 
range by Planck. We also note that for small enough m x , 
a CMB experiment allows us to probe the region of ther- 
mal cross-sections, and that Planck sensitivity will reach 
it, making it possible perhaps to find hints of particle 
DM in CMB data. 



We note that the constraints obtained here are several 
orders of magnitude more stringent than those obtained 
from the analysis of high-redshift proto-halos [HI (see also 
the recent 28]). We stress that our results apply also for 
standard models with s-wave annihilations, where the an- 
nihilation cross section does not depend on v. In this 
case, our results can be directly compared with the con- 
straints from astrophysical observations in the local uni- 
verse 1291. 
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V. CONCLUSIONS 

We have studied the effects of WIMP DM self- 
annihilation on recombination, looking for signatures in 
the CMB anisotropy and polarization. Our analysis has 
been performed under the assumption that the shower 
produced by the WIMP annihilation interacts "locally" 
with plasma, and a fraction / of the energy is absorbed 
on-the-spot by the baryons, contributing to its ioniza- 
tion. Our methodology is consistent with other analysis 
on the effect of decaying or self-annihilating, low-mass 
dark matter performed in the past. We examine a range 
of higher WIMP masses (1 GeV/c 2 -10 TeV/c 2 ), and find 
that current WMAP data already allow us to put in- 
teresting constraints on self-annihilation cross sections 
higher than the "standard" thermal value, in the range of 
those required to explain the PAMELA and ATIC data 
in terms of dark matter. Physically motivated by the 
very low relative velocity of DM particles at the time 
of recombination, our constraints on "Sommerfeld" en- 
hanced cross sections are competitive with the existing 
ones from local Universe (galactic substructures), and 
an independent test achieved with standard physics of 



the early Universe. By using simulated mock data, we 
have found that the expected enhanced sensitivity of the 
Planck mission will improve our capability to constrain 
Sommerfeld enhancement in dark matter particle models, 
thus hinting toward the exciting possibility to be finding 
traces of particle dark matter in future CMB data. In- 
terestingly, Planck will very likely be able to probe the 
region of the thermal annihilation cross-section for low 
WIMP masses (<50GeV/c 2 ), the actual value depending 
on the gas-shower coupling /. Ultimately, a cosmic vari- 
ance limited experiment permits the possibility to probe 
cross-sections values at the order of the thermal one for 
0(100 GeV/c 2 ) mass WIMPs. 
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